Process Time Distribution of Driven Polymer Transport 
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We discuss the temporal distribution of dynamic processes in driven polymer transport inherent 
to flexible chains due to stochastic tension propagation. The stochasticity originates from the 
disordered initial configuration of an equilibrium polymer coil, which results in random paths for 
tension propagation. We consider the process time for when translocation occurs across a fixed pore 
and when stretching occurs by pulling the chain end. A scaling argument for the mean and standard 
deviation of the process time is provided using the two-phase picture for stochastic propagation. 
The two cases are found to differ remarkably. The process time distribution of the translocation 
exhibits substantial spreading even in the long-chain limit, unlike that found for the dynamics of 
polymer stretching. In addition, the process time distribution in the driven translocation is shown 
to have a characteristic asymmetric shape. 

PACS numbers: 36.20.Ey,87.15.H-,83.50.-v 



I. INTRODUCTION 

Polymer transport is a ubiquitous and fundamental 
process in biological science and technology. A good ex- 
ample is polymer translocation in which a polymer (e.g., 
DNA or RNA) is driven across a small pore by a chem- 
ical or electric potential difference. Individual translo- 
cation events can be detected by applying an external 
force [H-01: allowing the fundamental properties of poly- 
mer transport to be investigated. In addition, many the- 
oretical investigations have been conducted [8 : 18]. Sta- 
tistical results obtained in experiments exhibit a wide 
distribution in translocation times (i.e., the time taken 
for a polymer to pass through a pore) 0, @. In addi- 
tion, numerical simulations have generated considerable 
insight [l3|, [HI, Il8rl26| , which also indicated a broad dis- 
tribution of the translocation time [25|, [2(| . Currently, 
the mean translocation time is attracting a lot of interest, 
particularly its scaling exponents with respect to chain 
length Nqcl (No is the segment number with size a) and 
driving force /. However, fluctuations in the transloca- 
tion time need to be clarified to fully understand the 
translocation process, but they have received less atten- 
tion to date. 

As is well known, a polymeric chain in equilibrium is 
a fluctuating coil whose configuration is mathematically 
described by the trajectory of a random (or self-avoiding) 
walk. It seems reasonable that when a force is suddenly 
applied to a portion of the chain, the entire chain will 
not immediately respond to the stimulus, but only the 
portion of the chain that is initially set in motion. The 
responding domain gradually evolves over time with the 
propagation of tension along the chain backbone. In do- 
main growth, the propagation front of the tension fol- 
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lows the disordered trajectory characteristics of a ran- 
dom coil configuration. The present study investigates 
the stochastic dynamics of tension propagation, which 
affects the process time distribution of driven polymer 
transport. 



II. ORIGIN OF STOCHASTICITY 

There are generally two sources of the stochasticity in 
stochastic processes: noise (random forces) and uncer- 
tainty in the initial distribution. In the case of polymer 
transport, the former source is the Brownian force, which 
is thermal in origin, while the latter source is reflected in 
the equilibrium configuration when one end of the chain 
arrives at the pore at time t = 0. To quantify the relative 
importance of these two sources, we consider the time 
scale involved in polymer transport over a distance X by 
an external force. Two time scales can be identified: the 
diffusion time to(X) ~ X 2 /D, where D is the diffusion 
coefficient, which depends on the chain conformation, 
and the convection time Tf(X) ~ X/V, where V = f/T 
is the mean biased velocity and T is the frictional coeffi- 
cient, which is related to the diffusion coefficient through 
the Einstein relation D = k B T/T. The ratio Pe(X) = 
t d (X)/t s (X) ~ fX/(k B T) defines a Peclet number. 
From the condition Pe(R cq ) > 1, where R eq is the equi- 
librium coil size, the condition for the driven transport 
regime can be derived as fR eq /(k B T) > 1, where the ef- 
fects of thermal fluctuations become relatively unimpor- 
tant. For stronger forces Pe(a) > 1 fa/(k-QT) > 1, 
the Brownian force can be neglected even when consider- 
ing the transport of a single segment. On the other hand, 
there is inevitably an uncertainty in the initial distribu- 
tion, irrespective of the driving force. This uncertainty 
is expected to be a major cause of the wide translocation 
time distributions observed in experiments. 

Below, we neglect minor effects associated with the 
stochasticity due to the Brownian force. We consider 
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FIG. 1: (Color online) Schematic representation of driven 
polymer translocation (in the trumpet [TP] regime). The 
tensed moving domain is shaded. The propagation front of 
the tension is located at x — —Rn- The top and bottom snap- 
shots depict stochastic evolution along the initial disordered 
configuration during time interval At. 



segment length a, force k^T/a, and time rja 3 / {k-^T) , re- 
spectively (where k&T is the thermal energy and 77 is the 
solvent viscosity). 

Here, we introduce the relevant exponents. The Flory 
exponent v and the dynamical exponent z are associated 
with the static and dynamic properties, respectively; the 
equilibrium coil size is described by i? eq ~ Nq and the 
hydrodynamic drag force acting on the coil with moving 
velocity V is described by ~ Rl^ 2 V [27[. The expo- 
nents have the following values in practical important 
cases: v = 1/2 (for an ideal chain), v = v 3 ~ 0.5876... 
(for a three-dimensional self- avoiding chain), z = 3 (for 
non-draining case), and z = (1 + 2v)/v (for free-draining 
case). 

The moving domain is characterized by the dynamical 
equations of state, [TJ EH, EE, US HH, which describe 
the global polymer conformation in terms of the driv- 
ing force /, the extension Rn along the x-axis, and the 
representative velocity V: 

N(t)-M(t)=a F (f)R N , (1) 



the transport time distribution due to stochastic tension 
propagation along a random backbone trajectory for two 
cases: driven translocation across a pore for which the ex- 
ternal force is spatially fixed and stretching by suddenly 
pulling one chain end, for which the forcing point is fixed. 
The stochastic scenario allows higher cumulants (e.g., the 
standard deviation and skewness) to be discussed in ad- 
dition to the average of the process time distribution. 



III. DRIVEN TRANSLOCATION 

We first consider driven translocation. To analyze the 
nonequilibrium response, we adopt the two-phase pic- 
ture [HI, Ell Ell, Gil * n w hi cn the chain separates on the cis 
side into quiescent and moving domains. As illustrated 
in Fig. 1, the coordinate x is taken to be perpendicular 
to the wall with the pore. The driving force acts only 
at the pore site [x = 0) with a constant force magnitude 
/ in the x-direction from the cis to the trans side. The 
linear polymer consisting of No segments is initially in 
x < 0. The polymer segments are numbered from one 
end of the chain to the other end (-/Voth segment); the 
first segment arrives at the pore at time t = 0. As time 
evolves, the rear segments are gradually sucked and the 
moving domain on the cis side grows in the negative x- 
rcgion. At time t, M(t) is the number of the segment 
at the pore (x = 0) and N(t) is the segment number at 
the tension-propagation front (x = —Rn). The propaga- 
tion front follows a random path sampled from equilib- 
rium coil configurations. Consequently, the dynamics of 
a particular sample strongly depends on how the poly- 
mer is brought to the pore. Below, the length, force, and 
time are made dimcnsionless by employing units of the 
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where crp(f) is the line segment density at the fore end 
(pore) and /j ~ A "^, /* ~ 1, and /** ~ Nq are the 
characteristic forces separating different regimes. The 
derivation of these forces and characteristic regimes with 
deformed shapes [TP], [SF], and [SS] are described in 
Ref. [l5|. The deformation characteristics illustrated in 
Fig. 2 are briefly introduced below: (a) In the trum- 
pet [TP] regime, the nonuniform deformation is ana- 
lyzed in terms of the space-dependent blob model; (b) 
In the stem-flower [SF] regime, the front end is almost 
fully stretched, whereas the rear end retains a blob-like 
shape; (c) All the tensed segments are almost completely 
stretched in the strong-stretching [SS] regime. Check- 
ing the Peclet number, we find that Pe(R cq ) > 1 <^> 
/ ~ Nq", which covers the deformed shape regimes. 



A. Stochastic Evolution of Tension Propagation 

In this section, we discuss the stochastic time evolu- 
tion of tension propagation. By putting the segment flux 
at pore dM/dt = <Jf{I)vf into the time derivative of 
eq. ([T]), we have 



dN dR N 
— — = a F v F + op — 
dt dt 



(5) 
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obtained from J Q N R k dR k = ££i(l/2)(ifc+ii*-i)[ifc- 
i? fc _i] = (1/2) [i? 2 . - i? 2 ,] and the initial condition i? (t = 
0) = 0. This leads to the mean value: 

(vfZ)- 1 (R k }dk-(2Z)- 1 (R 2 N ) (12) 
Jo 

A 1 (\=l)(a F Z)- 1 N 1+v - Z~ 1 N 2u , (13) 

where (• • •) denotes the ensemble average over the initial 
conditions. To derive cq. (|13p. we assume that the spatial 
distance from the pore site to the fcth segment's position 
at the initial time (t = 0) has a normalized distribution 
function SHI: 

P 1 (R k ) = k- u ip 1 (R k /k v ). (14) 

This gives the mean value of the positive power of R k as 
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FIG. 2: (Color online) Illustration of the deformed shapes 
associated with the different driving force magnitude, (a) 
Trumpet (/»</</*); (b) stem-flower (/* </</**); (c) 
strong-stretching (/** < /). (See Ref. [15| for more details.) 



where v F — v(x — 0) is the segment velocity at the pore. 
To proceed, we apply the steady-state approximation by 
setting V = v F ; 



dTV di?AT 



(6) 



Note that the adopted steady-state ansatz V = v F is dif- 
ferent from that in our previous studies V = «r = v(x = 
—Rn) [ULIU, Ell- The reason for the present choice lies 
in the fact that it is consistent with the iso-flux condi- 
tion required for the moving domain in the translocation 
process (l6l | , while the previous one is not (see Appendix 
A). 

Equation ([6]) can be arranged into a stochastic differ- 
ential equation that expresses the stochastic evolution of 
t as a function of TV 



dt = AdN + BdR N 

A = l/(a F V) 
B = -1/V 



(7) 

(8) 
(9) 



where di?jy is the noise originating from the initial ran- 
dom configuration. 

At first sight, eq. (O seems to consist of a deterministic 
drift term and a noise term. However, this is not the case, 
since V is also a stochastic variables. By substituting 
eq. ([2]) into eq. ([7]), we obtain 

dt = {R N /(a F Z))dN - {R N /Z)dR N (10) 

Integrating this equation with respect to TV yields 



t=(a F Z)- 1 R k dk-{2Zy 1 R 
Jo 
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[11) 



where we adopt a midpoint discretization scheme in the 
definition of the stochastic integral 1 . The second term is 



1 We adopt the Stratonovich (midpoint discretization) scheme 
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dR fc = Ai(A)fc A ", (15) 



where the numerical coefficient Ai(A) = u^ipi (u)du 
follows from the variable transformation u = R k jk v . The 
fact that ipi( u ) is expected to decrease exponentially at 
the large u guarantees the convergence of the above in- 
tegral for positive A. We then substitute cq. (|15[) into 
cq. (fT2")) because the exponents of all terms in eq. (fT2"|) 
are positive. For a strong force / > 1 or a long chain 
TV ^> 1, the first term in eq. (|13j) is dominant, so we 
obtain the following scaling for the mean translocation 
time 
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(17) 



(18) 



where, as before, we assume a normalized joint distribu- 
tion function P 2 (R k ,R k ') = k-"k'-^ 2 (R k /k v \R k , /k' v ) 



because it is expected to be suitable for describing real physi- 
cal phenomena (random coil configuration) in terms of idealized 
delta-correlated noise and because it is compatible with conven- 
tional formulas in real analysis l30j . Note, however, that the 
scaling exponents are not altered in this case, even if the Ito 
integral is employed. 
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for the spatial distances Rk and Rk' at the initial time 
(t = 0). This enables the mean value of the product of 



R% and R%, to be calculated: 



where the 

fOO 



numerical coefficient 
'uV v - 



/A v 



(19) 



A a (A,A') = 

J Q du J du'u x u' x ip2 (u, v!) follows from the vari- 
able transformation u = Rk/k v and v! = Rk'/k' v . 
Combining eqs. (fT3|) and (IT7l) gives the standard devi- 
ation SD(t) = J (t 2 ) — (t) 2 . Since the cross-correlation 
{AR X AR X } = [A 2 (A,A) - Ai(X) 2 ]k x "k' Xu > 0, where 
AR X = R^ — (R^), we obtain the following scaling 
relation for the standard deviation 



SD(t) 



N 1 - 



Mf)z(fV 



(20) 



which has the same scaling structure as that of the mean 
time; i.e., SD(i) - (t). 

There are two distinct stages in the driven transloca- 
tion process: tension propagation and post-propagation. 
The above-mentioned scenario considers the tension- 
propagation stage only. In the post-propagation stage, 
the moving domain consists of the entire chain on the 
cis side [ll [H [HI 
then given by r = , p 
ond terms on the right-hand side respectively correspond 
to the tension-propagation and post-propagation peri- 
ods. Equation (j6]) can be applied to the post-propagation 
stage (N = Nq = const.) dynamics by using the trivial 
condition dNo/dt = 0. Combining it with eq. (HJ) leads 
to 



The total translocation time is 
n + T pp , where the first and sec- 



dt 



-V^dR 



N„ 



-ZifyiR^dRNo, (21) 



where it is noted that the stochasticity in the tension 
propagation pathway is no longer relevant, thus, eq. ([2Tj) 
is a deterministic equation. Solving this, we find the 
post-propagation time (r pp ) = Z{f)~ 1 (Rn ( ,(t p ) 2 ) ~ 
Z(f)~ 1 No 2u . This means that, at the scaling level, t pp is 
a correction term, so that the translocation time is even- 
tually written as r ~ r p [111, [l2|, EH . In this study, we 
regard the translocation time as r ~ r p . 



B. Discussion 

We have argued that the scaling exponents of the mean 
time are the same as those of the standard deviation 
in their respective regimes. To test the above scaling 
predictions, we numerically integrated the stochastic dif- 
ferential equation (eq. <jTUJ) ) for a Rouse chain (i.e., no 
excluded volume v = 1/2 and no hydrodynamic in- 
teractions z = 4). In this simplest case, the appar- 
ent random force dR^/dN becomes white noise [32( as 
dN f ,' } — 6(N — N'), which enable us to treat the 
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FIG. 3: Numerical calculation of the stochastic translocation 
process under the SF or SS regimes. Double logarithmic plots 
of the mean and standard deviation against (a) the segment 
number N and (b) the magnitude of the driving force /. 
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FIG. 4: (a) Frequency distribution of the process (tension- 
propagation) time t p = t(No,f) (~ r) for translocation with 
No = 500, / = 100. (b) Probability density function as a 
function of normalized propagation time (r p — (r p ))/SD(r p ) 
for various conditions. The SF or SS regimes are adopted and 
the sample number is 10, 001. 



Wiener process [31j . Such noise is generated by a reflect- 
ing boundary (x = 0) that the propagation front cannot 
cross. Figures 3-5 show the corresponding numerical re- 
sults in the SF or SS regimes (eq. © and bottom lines of 
eqs. Q and Q). The triangles in Fig. 3 indicate the the- 
oretical exponents; they show that the slopes correspond 
very well with the numerical results. 

The distribution function contains more information 
about the translocation time. The profiles obtained in 
most numerical studies are asymmetric having a sin- 
gle left-skewed peak {24U261 ] . Rapid translocation ex- 
periments have also typically exhibited asymmetric pro- 
files [H, Q, although some have given right-skewed or 
symmetric profiles 0, 0] ■ Figure 4 (a) demonstrates a 
histogram of the translocation time for the chain length 
A*o = 500 obtained by our numerical integration. A char- 
acteristic profile with a single left-skewed peak in the 
time distribution is clearly observed. Similar asymme- 
try was observed over a broad range of parameters, as 
shown in Fig. 4 (b) (where the translocation time is nor- 
malized as (t— (i))/SD(i)). To gain a better understand- 
ing of the asymmetry, the following skewness (=third 
cumulant/SD 3 ) was analyzed: 



noise term in the stochastic differential equation as the 
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FIG. 5: Plots of the normalized standard deviation (NSD) 
SD(i)/(t) against (a) the segment number N and (b) the mag- 
nitude of the driving force / under the SF or SS regimes. 



(t 3 )-3(f 2 )(t)+2(t)- 
SD(t) 3 



(22) 



Introducing the normalized joint distribution 
for the three segments Pz(Rk, Rk' , Rk") = 
k- v k'- v k"- v i)^{u,u' ,u") with u = R k /k v , u' = R k '/k' v , 
and u" = Rk"/k" u , we obtain (t 3 ) ~ {t) 3 through a 
similar discussion of the derivation of (t) and (t 2 ). Thus, 
the skewness 71 is constant on the scaling level and 
the profiles approximately overlap each other. Recent 
numerical results [24l l25j are consistent with the present 
finding. 

What factor is responsible for the asymmetric profile? 
An equilibrium coil in free space retains a symmetric con- 
figuration. However, this is not the case for the initial 
configuration of a translocating polymer when a wall ex- 
cludes the polymer, one end of which is located at the 
pore. To clarify this point, we check the initial segment- 
distribution function obtained as follows: 



P ini (AT,x) = y^VA e -^ 2JV \ 



(23) 



where x < and j° oo Pi n i(N,x)dx = N. This is ob- 
tained by solving the diffusion equation dNPi nl (N, x) = 
(l/2)d 2 Pi D i(N, x) with the reflecting boundary condi- 
tion d x Pi n i(N, x)\ x= q = 0. In cq. (|2"3"|) . the reversal 
point x re{ ; i.e., P ini {N,x Icl + x) = P ini (N,x Icl - x) for 
x IC { < docs not hold. In other words, the configura- 
tion symmetry is violated; this is the only source of the 
asymmetric translocation time distribution in our model. 
Since the tension-propagation mechanism under a strong 
force has received strong support from numerical stud- 
ies [H, [H, [24[, and, according to Fyta et al. [23j, the 
translocation time depends on the initial configuration, 
it is reasonable to assume that stochastic evolution of the 
tension propagation along the initial asymmetric config- 
uration gives rise to the left-skewed profiles observed in 
numerical and experimental studies. Conversely, the de- 
tection of other translocation time profiles (e.g., right- 
skewed and symmetric profiles [H) may be regarded as a 
good indicator that other factors (such as specific inter- 
actions in the pore) dominate the process. 

Figures 5 (a) and (b) show plots of the normalized stan- 
dard deviation (NSD) SD (£)/(£). It is nearly constant 



against the parameters / and Nq{~^> 1). The experimen- 
tal results in ref. Q collecting simple translocation events 
without folding show that the NSD has values in the 
range ~ 0.1 — 0.2 independent of the chain length; this is 
in qualitative agreement with our prediction. 



IV. STRETCHING 

In this section, we consider a different protocol for 
chain transport in which one end of a linear polymer 
chain is pulled suddenly, as shown in Fig. 6. If the lin- 
ear polymer chain is initially relaxed in solution and one 
end is pulled by a constant force / in the x-direction af- 
ter t = 0, the chain will be progressively deformed and 
will eventually settle in a steady state [34j . Applying the 
two-phase picture and the stochastic method in a simi- 
lar manner as above, we analyze the stretching dynamics 
and derive the mean time of the transient process and 
its standard deviation. The segments are numbered from 
the pulled segment (one chain end) to the other end. The 
positions of the front and the rear ends in the moving do- 
main are x = l(t) and —R;sr(t) respectively; the pulled 
segment is initially located at the origin. Here again, we 
utilize the dynamical equation of state to characterize the 
global conformation [341 ] : 



N = a F (f)L 
LV = Z(f). 



(24) 
(25) 



where the extension along the x-dircction is L = I + Rn 
(see Fig. 6), and ctf(/), 2(f) are given in eqs. ©, (@|, 
respectively. Evidently, eqs. (|2"4")l and (|2"5)) exactly cor- 
resonds to eqs. (HJ, respectively. 



A. Stochastic Evolution of Tension Propagation 

Combining the boundary condition dl(t)/dt = vp{t) — 
V(t) with the time derivative of I = L — Rn, we obtain 



dt 



dR 



N 



dt 



(26) 



f 

N - — I — I 1 — ►* 

-R N I 

FIG. 6: (Color online) Schematic representation of the 
stretching process (in the trumpet [TP] regime). The tensed 
moving domain is shaded. 
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Combining cqs. ([24]) and (|25|) with cq. ([26]) gives the 
stochastic equation: 



di = „4d/V + £d.R 



N 



(27) 



where di?w acts as an apparent random force with 
(di?jy) = 0, and 



A = a F (/)" 2 Z(/)- 1 7V 



(28) 
(29) 



In contrast to eqs. © and © for the translocation dy- 
namics, eqs. (|28| and (|29[) contain only deterministic 
variables. Integrating eq. (f2~T)) with respect to the seg- 
ment number leads to 



t 



N' 2 



1 



2eriZ opZ 



N 



dfc 



(30) 



We then take the ensemble average to give 
N 2 



2a F {fYZ{f) 

N 2 fi- Z (/a </</*)••■ [TP] f ^ 
iV 2 /- 1 (/* </)••• [SF], [SS] -^ iJ 



The mean square of the process time is calculated to be 

= {2a 2 F Z)~ 2 N 4 + c{a F Z)- 2 N 2+2l/ , (32) 

where c is a numerical coefficient of order unity. Note 
that (^^) = 5(k-k>) and (^^f) * {k~k') 2 ^ 2 
for ideal and self-avoiding chains, respectively. We then 
obtain the standard deviation 



SD(i) 



N 1+v 



Mf)z{fY 

N i+, f i +i - z </</*)... [TP] 

N^f- 1 (/* </)••■ [SF], [SS] 



(33) 



and the NSD 



SD(t) 



a F (/)iV 



(34) 



"Stretching" vs. "Translocation" 



Comparison of cqs. (|T6|) and ((20)) with eqs. (|31j) 
and ([33]) reveals that the process distributions for 
"translocation and "stretching" differ remarkably. In 
the former case, the scaling exponents of the mean and 
the standard deviation are identical so that SD (£)/(<) — > 
const. > in the limit N — > oo. In the latter case, 
cq. (J3JJ) indicates SD(t)/(<) -> in the limit N -> oo. To 
determine the cause for this difference between the two 



cases, we compare the relevant equations by listing them 
again: 



Translocation 

N - M = a F {f)L 
LV = Z(f) 
(TV — M)V = ap{f)Z{f) 



<^ Stretching 

r N = a F (f)L 

& { LV = Z(f) (35) 
I NV = a F (f)Z(f) 



R 



N 



& L = R 



N 



(36) 



Equations (|35[) are a set of different forms in the dynam- 
ical equations of state, eqs. (|3o| give the chain extension 
in the moving domain. The scaling forms of cr F (f) and 
Z(f) are given in eqs. (J3J) and (U]). As is clear from this 
comparison, the essential difference is found in eq. (|36|) . 
i.e., in the translocation process, the segments with their 
label smaller than M(t), which are in the trans-side, are 
already free from the driving tension, thus, are not rele- 
vant to the moving domain. 

First, this leads to the so-called iso-flux condition 
N(t) ~ M(t) to leading order One can indeed 

check N(t) ^> a F (R]sr(t)) using the first moment rela- 
tion (R N {t)) ~ N{tf (eq. (USD). 

Second, this causes the qualitative difference in the 
stochastic effect in two processes. In the stretching dy- 
namics, the evolution of the moving domain has two com- 
ponents (Fig. 6, the right equation in eq. (|36[) ). One com- 
ponent is the pulled out part in the region < x < I and 
the other component is the rear part (—Rn < x < 0) 
due to the tension spreading. The former dominates the 
process (I ^> Rn) giving rise to the deterministic term in 
the stochastic differential equations (eqs. (|2"T)) and (|28p). 
while the latter acts as a small noise term (eqs. ([2T]) 
and ([29]) ). In contrast, in the translocation dynamics, 
the moving domain has only a stochastic component R 
(Fig. 1, the left equation in eq. ([35|) ~). which results in the 
non- vanishing NSD, even in the long-chain limit. 



V. REMARKS AND SUMMARY 

Before concluding, let us make some comments con- 
cerning the relation between the present study and other 
theoretical studies of polymer translocation. 

First, we have focused on the "driven regime" (i.e., 
/ > /jj ~ Nq") for which the chain deformation dynam- 
ics are important alongside tension propagation. In par- 
ticular, as noted in Sec. II, we expect that our approach 
can accurately calculate the process time distribution in 
the strong force regime Pe(a) > 1 «=>/>/* ~ 1. In such 
situations, the retardation effect due to segment accumu- 
lation on the trans side provides only a weak perturbation 
so that it can be neglected for asymptotic scaling. How- 
ever, it is essential to appropriately treat the trans side 
to describe the translocation process in the weak-force 
"near-equilibrium" re gim e / < /j, as has been done in 
several studies H [TH, Tlij . In this regime, the Brownian 
force also greatly influences the process time distribution. 
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It is interesting to determine the scaling structure of the 
standard deviation in such a weak-force regime. 

Second, the pore has been assumed to be only geo- 
metric constrictions. We have restricted ourselves to this 
simplest situation to get a clear-cut impact of the initial 
polymer configuration on the process time distribution. 
However, elucidating the case with functionalized pores 
with friction, confinement or specific interactions with 
polymers is an important problem. 

To summarize, we have discussed the process time dis- 
tribution for driven polymer transport inherent in a flex- 
ible molecular chain. The tension propagates with time 
along random paths following the initial configuration 
and the stochastic propagation mechanism is introduced 
as the distribution origin. We give the formulation for the 
two cases of translocation and stretching, which involve 
different forcing points. The forcing point of the former 
process is fixed in the rest frame, whereas that of the lat- 
ter process is fixed to the chain end. Our analysis predicts 
the scaling exponents of the first and second cumulants 
(mean, variance) in the process time. In translocation dy- 
namics, the probability distribution of the process time 
has a characteristic asymmetric shape, which may reflect 
the initial shapes, and fluctuation effects remain substan- 
tial even in the long-chain limit. On the other hand, the 
broadness of the distribution in stretching dynamics be- 
comes unnoticeable in the long-chain limit. 



concreteness, we focus on the TP regime and investigate 
the velocity difference inside the moving domain under 
iso-flux condition. The flux balance opvp — &rvr is 
transformed as 



v F /v R = a R /a F ~ (1 - 5v R /V)- q UR) q 



(37) 



where vp = V + 5vp, vr = V — 5vr and <7r ~ gn/^R 



with q = (1 — v)/[{z — l)v]. Then we find 



i + ^ = w 
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(38) 



Note that (fR) q > 1 holds under the deformed shape 
regimes. 

(i) If we assume V(t) = Vr, the iso-flux condition 
eq. ([38]) becomes Sv/V = (fR) q — 1, which might not 
satisfy the condition for the steady-state approximation. 

(ii) If V(t) = vp, the iso-flux condition eq. (|3"5]l be- 

_ q 

comes Sv/V = 1 — (fR) 1_ « < 1, thus the condition for 
the steady state approximation is satisfied. In a similar 
way, the validity of the steady-state approximation with 
V = vp is verified in the SF regime, too, while the iso-flux 
condition is trivially realized under the SS regime. 

In addition, we can check the correspondence at the 
level of the dynamical equation. The local force balance 
under TP regime is generally given by 



f 
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APPENDIX A: STEADY STATE ANSATZ AND 
ISO-FLUX CONDITION 

Recently, Rowghanian and Grosberg pointed out a 
unique property inherent to the moving domain in the 
translocation process, that is to say, the segment flux is 
almost constant everywhere in the moving domain fl6| . 
This leads to the balance of the fluxes into and out of the 
moving domain, i.e., opVp ~ cjrvr, where vr = v(x = 
—Rn) and or = <j(x = — Rn)- The subscript "i2" indi- 
cates the rear part of the moving domain, where the last 



tensed blob of size £r 



.V(i-*) 
J R 



is located. In the ap- 



pendix, we shall make a connection between the iso-flux 
condition and our steady-state ansatz. 

The steady state approximation would be valid un- 
der the condition; r rc i ax 7 < 1 Sv < V with the rep- 
resentative (or average) velocity V(t) = v(x,t) with 
3x € [0, — Rn], the relaxation time r re i a x — Rn/V and 
the shear rate 7 ~ Sv/R = (vp — vr)/Rn- For the 



v{x,t)£,(x,t) z z dx 



(39) 



Its spatial derivative leads to —d x £,(x,t) 
v(x, t)£(x, t) 2 ' 1 . We find its solution as 



v(x, t)dx 



--Ra 



(40) 



Putting the force balance at pore / = 1/£f 
0, t) leads to 



Rn 



I 

Rn 



v(x, t)dx 



z{f). 



(41) 



This is the dynamical equation of state (eqs. ([2])) with 
the "average velocity" V{t) = J_ R v(x,t)dx. 

In evaluating the average velocity V(t), the fore end 
has the dominant weight under the iso-flux J{t) condi- 
tion, since v(x, t) = J(t)/a(x,t) increases towards the 
fore part given the nonuniform spatial segment profile. 

Thus, V(t) = [^J° RN v(x,t)dx\ ~ v(x = 0,t) = 

J(t)/cr(x = 0, t) is expected so that the reasonable corre- 
spondence between the steady state approximation with 
the ansatz V(t) = vp(t) and the iso-flux model [l6| is 
found. By multiplying up to both sides of eq. ([2]), we 
obtain the dynamical equation of state in terms of the 
flux J(t) as RpfJ — ap(f)Z(f), which coincides with 
the result based on the iso-flux model Hal. 
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